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Abstract Market makers continuously set bid and ask quotes for the stocks 
they have under consideration. Hence they face a complex optimization prob- 
lem in which their return, based on the bid-ask spread they quote and the fre- 
quency at which they indeed provide liquidity, is challenged by the price risk 
they bear due to their inventory. In this paper, we consider a stochastic con- 
trol problem similar to the one introduced by Ho and StoU 17 and formalized 
mathematically by Avellaneda and Stoikov '3^. The market is modeled using 
a reference price St following a Brownian motion with standard deviation a, 
arrival rates of buy or sell liquidity-consuming orders depend on the distance 
to the reference price St and a market maker maximizes the expected utility of 
its P&L over a finite time horizon. We show that the Hamilton- Jacobi-Bellman 
equations associated to the stochastic optimal control problem can be trans- 
formed into a system of linear ordinary differential equations and we solve the 
market making problem under inventory constraints. We also shed light on the 
asymptotic behavior of the optimal quotes and propose closed-form approxi- 
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mations based on a spectral characterization of the optimal quotes. 

Keywords Stochastic optimal control • High-frequency Market Making • 
Avellaneda-Stoikov problem 

1 Introduction 

From a quantitative viewpoint, market microstructure is a sequence of auc- 
tion games between market participants. It implements the balance between 
supply and demand, forming an equilibrium traded price to be used as refer- 
ence for valuation. The rule of each auction game (fixing auction, continuous 
auction, ...) is fixed by the firm operating each trading venue. Nevertheless, 
most of all trading mechanisms on electronic markets rely on market partic- 
ipants sending orders to a "queuing system" where their open interests are 
consolidated as "liquidity provision" or form transactions [T|. The efficiency 
of such a process relies on an adequate timing between buyers and sellers, to 
avoid too many non-informative oscillations of the transaction price (for more 
details and modeling, see for example |20|). 

In practice, it is possible to provide liquidity to an impatient buyer (re- 
spectively seller) and maintain an inventory until the arrival of the next im- 
patient seller (respectively buyer) . Market participants focused on this kind of 
liquidity-providing activity are called "market makers" . On one hand they are 
buying at the bid price and selling at the ask price they choose, making money 
out of this "bid-ask spread" . On the other hand, their inventory is exposed to 
price fiuctuations mainly driven by the volatility of the market (see [2.,5iil01 

The recent evolution of both technology and market regulation reshaped 
the nature of the interactions between market participants during continuous 
electronic auctions, one consequence being the emergence of "high-frequency 
market makers" who are said to be part of 70% of the electronic trades in 
the US (40% in the EU and 35% in Japan) and have a massively passive 
{i.e. liquidity-providing) behavior - a typical balance between passive and 
aggressive orders for such market participants being around 80% of passive 
interactions (see [12] )• 

From a mathematical modeling point of view, the market making problem 
corresponds to the choice of optimal quotes (i.e. the bid and ask prices) that 
such agents provide to other market participants, taking into account their in- 
ventory limits and their risk constraints often represented by a utility function 
(see [illllTllIIllMlllS]). 

Avellaneda and Stoikov proposed, in a widely cited paper [3 , an innovative 
framework for "market making in an order book" . In their approach, rooted to 
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an old paper by Ho and Stoll [T7] , the market is modeled using a reference price 
or fair price St following a Brownian motion with standard deviation a, and 
the arrival of a buy or sell liquidity-consuming order at a distance S from the 
reference price St is described by a point process with intensity Aexp{—kS), 
A and k being two positive constants which characterize statistically the liq- 
uidity of the stock. 

We consider the same model as in [3] - adding inventory limits - and 
we show, using a new change of variables, that the Hamilton- Jacobi-Bellman 
equations associated to the problem boil down to a system of linear ordinary 
differential equations. This new change of variables (i) simplifies the compu- 
tation of a solution since numerical approximation of partial differential equa- 
tions is now unnecessary, and (ii) allows to study the asymptotic behavior 
of the optimal quotes. In addition to these two contributions, we use results 
from spectral analysis to provide an approximation of the optimal quotes in 
closed-form. Finally, we provide in the case of our model with inventory limits 
a verification theorem that was absent from the original paper (the admissi- 
bility of the quotes obtained in the original Avellaneda-Stoikov model appears 
in fact to be an open problem!). 

Since Avellaneda and Stoikov seminal paper, other authors have consid- 
ered related market making models. Cartea, Jaimungal and Ricci [51 consider a 
more sophisticated model inspired from the Avellaneda-Stoikov onqj, including 
richer dynamics of market orders, impact on the limit order book, adverse se- 
lection effects and predictable a. They obtained closed-form approximations of 
the optimal quotes using a first-order Taylor expansion. Cartea and Jaimungal 
[7] recently used a similar model to introduce risk measures for high-frequency 
trading. Earlier, they used a model inspired from Avellaneda-Stoikov [5] in 
which the mid-price is modeled by a Hidden Markov Model. Guilbaud and 
Pham [14j also used a model inspired from the Avellaneda-Stoikov framework 
but including market orders and limit orders at best (and next to best) bid 
and ask together with stochastic spreads. Very recently, Guilbaud and Pham 
[15| used another model inspired from the Avellanada- Stoikov one in a pro- 
rata microstructure. 

It is also noteworthy that the model we use to find the orders a market 
maker should optimally send to the market has been used in a totally different 
domain of algorithmic trading: optimal execution. Bayraktar and Ludkovski 
P and Gueant, Lehalle, Fernandez- Tapia 12 used indeed a similar model to 
optimally liquidate a portfolio. 

This paper starts in section 2 with the description of the model. Section 3 is 
dedicated to the introduction of our change of variables and solves the control 
problem in the presence of inventory limits. Section 4 focuses on the asymptotic 
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behavior of the optimal quotes and characterizes the asymptotic value using 
an eigenvalue problem that allows to propose a rather good approximation 
in closed- form. Section 5 generalizes the model in two different directions: (i) 
the introduction of a drift in the price dynamics and (ii) the introduction 
of market impact that may also be regarded as adverse selection. Section 6 
carries out the comparative statics. Section 7 provides backtests of the model. 
Adaptations of our results are in use at Cheuvreux. 

2 Setup of the model 

Let us fix a probability space { f2, P) equipped with a filtration {Tt)t>o satis- 
fying the usual conditions. We assume that all random variables and stochastic 
processes are defined on {fl^T, (J"t)t>o,P). 

We consider a high-frequency market maker operating on a sing le stoclfl. 
We suppose that the mid-price of this stock or more generally a reference 
priceH of the stock moves as an arithmetic Brownian motiorQ: 

dSt = adWt 

The market maker under consideration will continuously propose bid and 
ask prices denoted respectively and S'f and will hence buy and sell shares 
according to the rate of arrival of market orders at the quoted prices. His 
inventory q, that is the (signed) quantity of shares he holds, is given by 

qt = - K 

where iV'' and N°- are the point processes (independent of (Wt)t) giving the 
number of shares the market maker respectively bought and sold (we assume 
that transactions are of constant size, scale^l to 1). Arrival rates obviously 
depend on the prices and quoted by the market maker and we assume, 
in accordance with the model proposed by Avellaneda and Stoikov [SJ, that 
intensities A** and A" associated respectively to and N'^ depend on the 
difference between the quoted prices and the reference price {i.e. 5\ = St — S\ 
and 6^ = — St) and are of the following forrrj^: 

A*(^'') = Ae-''^' = Aexp{-k{s-s^)) X''{S'') = Ae'''^'' = Aexp(-fc(s'^-s)) 

^ We suppose that this high-frequency market maker does not "make" the price in the 
sense that he has no market power. In other words, we assume that the size of his orders 
is small enough to consider price dynamics exogenous. Market impact will bo introduced in 
section 5. 

This reference price for the stock can be thought of as a smoothed mid-price for instance. 

* Since we will only consider short horizon problems, this assumption is almost equivalent 
to the usual Black-Scholes one. 

^ The only important hypothesis is the constant size of orders since we can easily replace 
1 by any positive size A. 

^ Some authors also used a linear form for the intensity functions - see 1171 for instance. 
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where A and k are positive constants that characterize the hquidity of the 
stock. In particular, this specification means - for positive 5^ and (5° - that 
the closer to the reference price an order is posted, the faster it will be executed. 

As a consequence of his trades, the market maker has an amount of cash 
evolving according to the following dynamics: 

dXt = {St + 5-i)dNf - {St - 5\)dN^t 

To this original setting introduced by Avellaneda and Stoikov (itself fol- 
lowing partially Ho and Stoll [H]), we add a bound Q to the inventory that 
a market maker is authorized to have. In other words, we assume that a mar- 
ket maker with inventory Q {Q > Q depending in practice on risk limits) will 
never set a bid quote and symmetrically that a market maker with inventory 
— Q, that is a short position of Q shares in the stock under consideration, will 
never set an ask quote. This realistic restriction may be read as a risk limit 
and allows to solve rigorously the problem. 

Now, coming to the objective function, the market maker has a time horizon 
T and his goal is to optimize the expected utility of his P&L at time T . In 
line with [3] , we will focus on C ARA utility functions and we suppose that the 
market maker optimizes: 

sup E [— exp (— 7(Xt -f- (/T'S't))] 

where A is the set of predictable processes bounded from below, 7 is the 
absolute risk aversion coefficient characterizing the market maker, Xt is the 
amount of cash at time T and qtSt is the evaluation of the (signed) remaining 
quantity of shares in the inventory at time T (liquidation at the reference price 



3 Characterization of the optimal quotes 

The optimization problem set up in the preceding section can be solved using 
the classical tools of stochastic optimal control. The first step of our reason- 
ing is therefore to introduce the Hamilton- Jacobi-Bellman (HJB) equation 
associated to the problem. More exactly, we introduce a system of Hamilton- 
Jacobi-Bellman partial differential equations which consists of the following 
equations indexed by q e {—Q, . . . , Q} for {t, s, x) e [0, T] x R^: 



^ Our results would be mutatis mutandis the same if we added a penalization term 
—b{\qT\) for the shares remaining at time T. The rationale underlying this point is that 
price risk prevents the trader from having important exposure to the stock. Hence, qt should 
naturally mean-revert around 0. 
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For \q\ < Q: 



dtu{t, X, q, s) + -a^dl^uit, x, q, s) 

+ sup X^{6^) [u{t, x — s + 6^,q + l,s) — u{t, x, q, s 



+ supA''(5") [u{t,x + s + 



(5", q — l,s) — u{t, X, q, s)] = 



For q 




+ supX''{S'')[u{t,x + s + S'',Q - l,s) - u{t,x,Q,s)] = 



For q 




+ supX''{6'') [u{t,x- s + 6\-Q + l,s)-u{t,x,-Q,s)] =0 



with the final condition: 



e {-Q, 



Q}, 



u{T, X, s) = - exp (-7(3; + g.s)) 



To solve these equations we will use a change of variables based on two 
different ideas. First, the choice of a CARA utility function allows to factor 
out the Mark-to-Market value of the portfolio {x + qs). Then, the exponential 
decay for the intensity functions A** and A" allows to reduce the Hamilton- 
Jacobi-Bellman (HJB) equations associated to our control problem to a linear 
system of ordinary differential equations: 

Proposition 1 (Change of variables for (HJB)) Let us consider a fam- 
ily {vq)\q\<Q of positive functions solution of: 

Vg e {-Q + 1, . . . , Q - 1}, Vgit) = aq\{t) - ri (ug-i(i) + Vg+i{t)) 

i>Q{t) = aQ'^vqit) - r/UQ_i(i) 

v-qit) = aQ'^v^qit) - T]v^Q+i{t) 

with Vg G {-Q, . . . ,Q},Vg{T) ^ I, where a = ^ja^ and ri = A(l + • 

_ 2. 

Then, u{t, x, q,s) = — exp(— 7(0; + qs))vq{t) is solution of (HJB). 

Then, the following proposition proves that there exists such a family of 
positive functions: 
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Proposition 2 (Solution of the ordinetry differential equations) Let us 

introduce the matrix M defined by: 



(aQ"^ -T] 

-r\ a(Q- 1)2 -ry ■•• 




M = 



\ 





-r;a(Q-l)2 -Tj 
••• -T] aQ^J 



V 

where a = ficr^ and rj = A{1 + ^)~^^^^''. 
Let us define 

V{t) = {V-Q{t),V-Q+-i_{t), . . .,Vo{t), UQ-l(i), VQ{t)y 

= eM-M{T-t)) X (1,...,1)' 

Then, {vq)\q\<Q is a family of 'positive functions solution of the equations of 
Proposition 1. 



Using the above change of variables and a verification approach, we are 
now able to solve the stochastic control problem, that is to find the value 
function of the problem and the optimal quotes: 

Theorem 1 (Solution of the control problem) Let consider {vq)\q\<.Q as 
in Proposition 2. 

Then u{t, x, q,s) = — exp(— 7(0; + qs))vq{t)~^ is the value function of the 
control problem. 

Moreover, the optimal quotes are given by: 

s - s- it, q, s) = 5"* it, q) = lln (^) + ^ In (l + J) , q^Q 

s-(i,g,s)-s = <5-(t,g)==^ln(^-^^^ +iln(l + |), q ^ -Q 
and the resulting bid-ask spread quoted by the market maker is given by: 
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4 Asymptotic behavior and approximation of the optimal quotes 

To exemplify our findings and in order to motivate tlie asymptotic approxi- 
mations we shall provide, we plotted on Figure [T] and Figure [5] the behavior 
with time and inventory of the optimal quotes. The resulting bid-ask spread 
quoted by the market maker is plotted on Figure |21 

We clearly see that the optimal quotes are almost independent of t, as 
soon as t is far from the terminal time T. This observation is at odds with the 
approximations propose^l in Avellaneda and Stoikov [3] using an expansion 
in q. It motivates however the study of the asymptotic behavior of the quotes. 




* In [3], the approximations obtained by the authors using an expansion in q leads to the 
following expressions: 

.-.iln(l-,2) + i±^,..(^_,) 

One can easily show, using the results of Theorem 1, that these approximations are nothing 
but the Taylor expansions of the optimal quotes for t close to T. 
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Inventory 

Fig. 2 Behavior of the optimal ask quotes with time and inventory, a = 0.3 Tick • s~^/^, 
A = 0.9 s-\fc = 0.3 Tick-i, 7 = 0.01 Tick-i,T = 600 s. 




600 



Inventory 



Fig. 3 Behavior of the resulting bid-ask spread with time and inventory, cr = 0.3 Tick • 
s-i/2, A = 0.9 s-i,fc = 0.3 Tick-\7 = 0.01 Tick"!, T = 600 s. 
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Theorem 2 (Asymptotics for the optimal quotes) The optimal quotes 
have asymptotic limits 

lim S''*{0,q)=6'^{q) 
T— S-+00 

lim <5-(0,Q)=C(g) 

T— i.+oo 

that can be expressed as: 



= i (l + 2) + 1 1„ (£] (,) = i In (l + 2) + 1 1„ (1^' 



where f'^ G M^'3+-'^ is an eigenvector corresponding to the smallest eigenvalue 
of the matrix M introduced in Proposition 2 and characterized (up to a mul- 
tiplicative constant) by: 

Q Q-i 
f e argmin ^ ag^/ + V J2 (Z^+i " + '^/q' + ^/-q' 

The resulting bid-ask spread quoted by the market maker is asymptotically: 

The above result, along with the example of Figure [TJ Figure [5] and Figure 
ini encourages to approximate the optimal quotes and the resulting bid-ask 
spread by their asymptotic value. These asymptotic values depend on /° and 
we shall provide a closed- form approximation for 

The above characterization of corresponds to an eigenvalue problem in 
]^2Q-i-i propose to replace it by a similar eigenvalue problem in 

for which a closed-form solution can be computed. More precisely we replace 
the criterion 



Q Q-i 
/" e argmin ^ ag^/g' + V " ^9)' + ^^q' + 



by the following criterion for G L^(M): 



+00 



/° e argmin / (^ax'^ f{x)^ + rjf (x)^^ dx 



The introduction of this new criterion is rooted to the following proposition 
that states (up to its sign) the expression for in closed form: 
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Proposition 3 Let us consider 

/° e argmin / (ax'^f{x)'^+r]f{x)'^jdx 

Then: 

From the above proposition, we expect to beliave, up to a multiplicative 
constant, as exp {^~\'J^<f^ ■ This heuristic viewpoint induces an approxima- 
tion of the optimal quotes and the resulting optimal bid-ask-spread: 



C(5) = il..(l + J)-^yf(2,-l) 




Fig. 4 Asymptotic behavior of optimal bid quote (bold line). Approximation (dotted line). 
Left: cr = 0.4 Tick ■ s-i/2_ A = 0.9 s-\ fc = 0.3 Tick"!, 7 = 0.01 Tick-\ T = 
600 s. Right: cr = 1.0 Tick ■ s-^/'^, A = 0.2 s-\ k = 0.3 Tick"!, 7 = 0.01 Tick"!, 
T = 600 s. 



We exhibit on Figure|l]and Figure 5 the values of the optimal quotes, along 
with their associated approximations. Empirically, these approximations for 
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Fig. 5 Asymptotic behavior of optimal ask quote (bold line). Approximation (dotted line). 
Left: a = 0.4 Tick ■ g-^/^, A = 0.9 s-\ k = 0.3 Tick"\ 7 = 0.01 Tick"!, T = 
600 s. Right: a = 1.0 Tick • s'^/'^, A = 0.2 s'^ , k = 0.3 Tick"!, 7 = 0.01 Tick"!, 
T = 600 s. 

the quotes are satisfactory in most cases and are always very good for small 
values of the inventory q. In fact, even though /° appears to be well approxi- 
mated by the gaussian approximation, we cannot expect a very good fit for the 
quotes when q is large because we are approximating expressions that depend 
on ratios of the form or -r^S— . 

Jg+l Jq-1 



5 Extensions of the model 

5.1 The case of a trend in the price dynamics 

So far, the reference price was supposed to be a Brownian motion. In what 
follows we extend the model to the case of a trend in the price dynamics: 

dSt = iJ,dt + adWt 

In that case we have the following proposition (the proof is not repeated) : 

Proposition 4 (Solution with a drift) Let us consider a family of func- 
tions {vq)iqi<:Q solution of the linear system of ODEs that follows: 

Vg e {-Q + 1, . . . , Q - 1}, vq{t) = {aq^ - (3q)vg{t) - r? {vg-i{t) + Vg+i{t)) 

VQ{t) - {aQ^ - l3Q)vQ{t) - rjVQ-iit) 
v-Q{t) = (aQ2 + pQ)v_Q{t) - v^-Q+iit) 
with Vg e {—Q, . . . ,Q},Vq{T) = 1, where a = ^^(7^, (3 = kjjL and rj = 
A(l + i)-(^+^). 

Then, u{t, x, q,s) = — exp(— 7(0; + qs))vq{t) is the value function of the 
control problem. 
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The optimal quotes are given by: 

s-s'^t,q,s) = 5'*{t,q) = hn(^^ 

k \Vq+l{t) 

q,s)-s = S"*{t, q) = Unf ""'^^^ 



k \Vg-l{t)J 7 

and the resulting bid-ask spread of the market maker is : 



r{t,q) 



k 



Moreover, 



^+00 ''7V kJ k 



lim 



lim (5"*(0,g) 



HI 



lim V'*(0, (?) = — T In 



1 , I fq+lfq-l 



■la 



where /° is an eigenvector corresponding to the smallest eigenvalue of: 
/aQ2 _ 



-V 







a(Q - 1)2 - /3(Q - 1) -r, 



V 







-rj a{Q - If - P{Q - 1) 







-V 



In addition to this theoretical result, we can consider an approximation 

similar to the approximation used for the initial model with no drift. We then 
obtain the following approximations for the optimal quotes and the bid-ask 
spread: 



C(.):^^ln(l + J) 



7Cr- 

JJ_ 
7(T^ 



M , 2g+l 
2 



2 

2q-l 



2kA \ k 
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5.2 The case of market impact 

Another extension of the model consists in introducing market impact. The 
simplest way to proceed is to consider the following dynamics for the price: 



When a limit order on the bid side is filled, the reference price decreases. On 

the contrary, when a limit order on the ask side is filled, the reference price in- 
creases. This is in line with the classical modeling of market impact for market 
orders, ^ being a constant since the limit orders posted by the market maker 
are all supposed to be of the same size. 

Adverse selection is another way to interpret the interaction we consider be- 
tween the price process and the point processes modeling execution: trades on 
the bid side arc often followed by a price decrease and, conversely, trades on 
the ask side are often followed by a price increase. 

In this modified framework, the problem can be solved using a change of 
variables that is slightly more involved than the one presented above but the 
method is exactly the same and we have the following proposition (the proof 
is not repeated): 

Proposition 5 (Solution with market impact) Let us consider a family 
of funetions {vq)\q\<Q solution of the linear system of ODEs that follows: 



e {-Q + 1, . . . ,Q - l},Vg{t) = aq^Vgit) - rje'^^ {vq-i{t) + Vq+i{t)) 



Then, u{t, x, q,s) = — exp(— 7(3; + qs + ^iq^))vq{t) *= is the value function 
of the control problem. 

The optimal quotes are given by: 



dSt = adWt + ^dN^ - ^dN^, 



^>0 




and T] = 





and the resulting bid- ask spread of the market maker is : 
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Moreover, 



lim S''*{0,q) = -ln(l + l)+i + Un(^] 
-^+00 7 V kJ 2 k \f^+ij 

lim 5-(0,g) = iln(l + J)+^ + iln(-^l 



lim V*(0,g) , 

T->-+oo k 



-In' 



w/iere /"^ is an eigenvector corresponding to the smallest eigenvalue of: 
( -?,e-*« \ 



-'qe 




2 ? 



V 



-r?e-4f a(g- 1)2 



-rye 



2 4 







—rje 



In addition to this theoretical result, we can consider an approximation 
similar to the approximation used for the initial model. Wc then obtain the 
following approximations for the optimal quotes and the bid-ask spread: 



<5^(g):..iln(l + 
C(a):^^ln(l + J)+f- 



C , 2(7+1 

— e*^ 

2 2 



2kA \ k) 



? 2g-l .J 



(727 
2fc]4 



1+^ 



■In 



6 Comparative statics 

Wc argued in section 4 that the value of the optimal quotes was almost inde- 
pendent of t for t sufficiently far from the terminal time T and we characterized, 
using spectral arguments, the asymptotic value of the optimal quotes. In or- 
der to obtain closed-form formulae, wc also provided approximations for the 
asymptotic value of the optimal quotes. Although these closed-form formu- 
lae are only approximations, they provide a rather good intuition about the 
influence of the different parameters on the optimal quotes. 
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6.1 Dependence on cr^ 

The dependence of optimal quotes on tr^ depends on the sign of the inventory. 
More precisely, we observe numerically, in accordance with the approximations, 
that: 



0cr2 


<o, 


9<5" 


>o, 


if 9 < 


da/' 


>o. 




>o, 


if g = 


>o. 


<o, 


if 9 > 



For the bid-ask spread, we obtain: 

The rationale behind this is that an increase in cr^ increases inventory risk. 
Hence, to reduce this risk, a market maker that has a long position will try 
to reduce his exposure and hence ask less for his stocks (to get rid of some 
of them) and accept to buy at a lower price (to avoid buying new stocks). 
Similarly, a market maker with a short position tries to buy stocks, and hence 
increases its bid quote, while avoiding short selling new stocks, and he increases 
its ask quote to that purpose. Overall, due to the increase in price risk, the 
bid-ask spread widens as it is well instanced in the case of a market maker with 
a flat position (this one wants indeed to earn more per trade to compensate 
the increase in inventory risk). 



6.2 Dependence on 

The dependence of optimal quotes on the drift ^ is straightforward and corre- 
sponds to the intuition. If the agent expect the price to increase (resp. decrease) 
he will post orders with higher (resp. lower) prices. Hence we have: 

ggb* gga* 

°° < °° > 



6.3 Dependence on A 



Because of the form of the system of equations that defines v, the dependence 
on A must be the exact opposite of the dependence on cr^: 
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For the bid-ask spread, we obtain: 



< 



dA 

The rationale behind these resuhs is that an increase in A reduces the 
inventory risk. An increase in A indeed increases the frequency of trades and 
hence reduces the risk of being stuck with a large inventory (in absolute value) . 
For this reason, an increase in A should have the same effect as a decrease in 

^2 



6.4 Dependence on 7 

Using the closed-form approximations, we see that the dependence on 7 is 
ambiguous. The market maker faces indeed two different risks that contribute 
to inventory risk: (i) trades occur at random times and (ii) the reference price 
is stochastic. But if risk aversion increases, the market maker will mitigate 
the two risks: (i) he may set his quotes closer to one another to reduce the 
randomness in execution and (ii) he may widen his spread to reduce price 
risk. The tension between these two roles played by 7 explains the different 
behaviors we may observe, as on Figure [S] and Figure [7] for the bid-ask spread 
resulting from the asymptotic optimal quotes: 
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6.5 Dependence on k 

From the closed-form approximations, wc expect to be decreasing in k for 
q greater than some negative threshold. Below this threshold, we expect it to 
be increasing. Similarly we expect 6^ to be decreasing in k for q smaller than 
some positive threshold. Above this threshold we expect it to be increasing. 



Eventually, as far as the bid-ask spread is concerned, the closed-form ap- 
proximations indicate that the resulting bid-ask spread should be a decreasing 
function of k. 

In fact several effects are in interaction. On one hand, there is a "no- volatility" 
effect that is completely orthogonal to any reasoning on the inventory risk: 
when k increases, in a situation where and S"" arc positive, trades occur 
closer to the reference price St- For this reason, and in absence of inventory 
risk, the optimal bid-ask spread has to shrink. However, an increase in k also 
affects the inventory risk since it decreases the probability to be cxeciited (for 
(5'', (5" > 0). Hence, an increase in k is also, in some aspects, similar to a de- 
crease in A. These two effects explain the expected behavior. 
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Fig. 8 Asymptotic optimal bid quotes for different inventories and different values of k. 
(7 = 0.3 Tick - ^ = 0.9 s"!, 7 = 0.01 Tick-\ T = 600 s. 

Numerically, we observed that the "no-volatility" effect dominated for the 
values of the inventory under consideration (see Figure \E\ for the case of the 
bid quot^. 

6.6 Dependence on the market impact ^ 

The market impact introduced in 5.2 has two effects on the optimal quotes. In 
the absence of price risk, given the functional form of the execution intensities, 
the direct effect of ^ is approximately to add | the each optimal quote: the 
market maker approximately maintains his profit per round trip on the market 
but the probability of occurrence of a trade is reduced. This adverse selection 
effect has a side-effect linked to inventory risk: since adverse selection gives 
the market maker an incentive to post orders deeper in the book, it increases 
the risk of being stuck with a large inventory for a market maker holding 
such an inventory. As a consequence, for a trader holding a positive (resp. 
negative) inventory, there is a second effect inciting to buy and sell at lower 
(resp. higher) prices. These two effects are clearly highlighted by the closed- 
form approximations exhibited in the previous section: 





side— ef feet 



^ The case of the ask quote is obviously similar. 
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adverse selection 



7 Backtests 

Before using the above model on historical data, we need to discuss some fea- 
tures of the model that need to be adapted before any backtest is possible. 

First of all, the model is continuous in both time and space while the real 
control problem is intrinsically discrete in space, because of the tick size, and 
in time, because orders have a certain priority and changing position too often 
reduces the actual chance to be reached by a market order. Hence, the model 
has to be reinterpreted in a discrete way. In terms of prices, quotes must not 
be between two ticks and we decided to round the optimal quotes to the near- 
est tick. In terms of time, an order of size AT is sent to the market and is 
not canceled nor modified for a given period of time At, unless a trade occurs 
and, though perhaps partially, fills the order. Now, when a trade occurs and 
changes the inventory or when an order stayed in the order book for longer 
than At, then the optimal quote is updated and, if necessary, a new order is 
inserted. 

Concerning the parameters, cr, A and k can be calibrated on trade-by-trade 
limit order book data while 7 has to be chosen. However, it is well known by 
practitioners that A and k have to depend at least on the actual market bid- 
ask spread. Since we do not explicitly take into account the underlying market, 
there is no market bid-ask spread in the model. For the backtest example we 
present below, A and k have been chosen independent of the spread but, in 
practice, the value of A and k are function of the market bid-ask spread. As 
far as 7 is concerned, we decided in our backtests to assign 7 an arbitrary 
value for which the inventory stayed between -10 and 10 during the day we 
considered (the unit being the ATS). 

Turning to the backtests, they were carried out with trade-by-trade data 
and we assumed that our orders were entirely filled when a trade occurred at 
or above the ask price quoted by the agent. Our goal here is just to exemplify 
the use of the modeO and we considered the case of the French stock France 
Telecom on March 15*'' 2012. 



ATS is the average trade size. 

We, voluntarily, do not give full details about the algorithm based on the model. 
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We first plot the price of tlie stock France Telecom on March 15*'' 2012 
on Figure [SI the evolution of the inventor on Figure HU] and the associated 
P&L (the stocks in the portfolio are evaluated at mid-price) on Figure [TTJ 




15:00 



16:00 



Fig. 9 Price of the stock France Telecom on 15/03/2012, from 10:00 to 16:00. 
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1400 




10:00 11:00 12:00 13:00 14:00 15:00 16:00 

Time 

Fig. 11 P&L when the strategy is used on France Telecom (15/03/2012) from 10:00 to 
16:00. 

This P&L can be compared to the P&L of a naive trader (Figure [12]) who 
only posts orders at the first hmit of the book on each side, whenever he is 
asked to post orders - that is when one of his orders has been executed or 
after a period of time At with no execution. 



800 




-600' ' ' ' ' ' ' 

10:00 11:00 12:00 13:00 14:00 15:00 16:00 

Time 



Fig. 12 P&L of a naive market maker on France Telecom (15/03/2012) from 10:00 to 16:00. 

Now, to better understand the details of the strategy, we focused on a 
subperiod of 1 hour and we plotted the state of the market along with the 
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quotes of the market maker (Figure I13p . Trades occurrences involving the 
market maker are signahed by a dot. 




12:00 12:10 12:20 12:30 12:40 12:50 13:00 



Fig. 13 Details for the quotes and trades when the strategy is used on Prance Telecom 
(15/03/2012). Thin lines represent the market while bold lines represent the quotes of the 
market maker. Dotted lines are associated to the bid side while plain lines are associated to 
the ask side. Black points represent trades in which the market maker is involved. 



Conclusion 

In this paper we present a model for the optimal quotes of a market maker. 
Starting from a model in line with Avellaneda and Stoikov [3, and rooted to 
Ho and StoU p^, we introduce a change of variables that allows to trans- 
form the HJB equation into a system of linear ordinary differential equations. 
This transformation allows to find the optimal quotes and to characterize 
their asymptotic behavior. Closed-form approximations are also obtained us- 
ing spectral analysis. 

The change of variables introduced in this paper can also be used to solve 
the initial equations of Avellaneda and Stoikov [5] and we provide a complete 
mathematical proof in |13j . However, in the absence of inventory limits, no 
proof of optimality is available for the quotes claimed to be optimal in [3] and 
their admissibility appears to be an open problem. 

An important topic for future research consists in generalizing the model to 
any intensity function. This is particularly important because the exponential 
form of the intensity is only suited to liquid stocks with a small bid-ask spread. 
Another important topic consists in introducing "passive market impact" {i.e. 
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the perturbations of the price formation process by hquidity provision). This 
is a real modehng challenge since no quantitative model for this type of impact 
has been proposed in the literature. 
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Appendix: Proofs of the results 

Proof of Proposition 1, Proposition 2 and Theorem 1: 

Let us consider a family {vq)\q\<Q of positive functions solution of the 
system of ODEs introduced in Proposition 1 and let us define u{t, x, q, s) = 
- exp {-j{x + qs)) Vq{t)~^ . 

Then: 

Now, concerning the hamiltonian parts, we have for the bid part {q ^ Q): 
sup A^((5^) [u{t, X — s + S^,q + l,s) — u{t, x, q, s)] 

5<' 



= supAe u{t,x,q,s) 

The first order condition of this problem corresponds to a maximum (be- 
cause u is negative) and writes: 

(t + .)exp(-..'-)(5teia)"*=t 

Hence: 




and 

sup X^{S^) [u{t, X — s + S^,q + l,s) — u{t, x, q, s)] 

S'' 

= - /^ Aexp{-k5^*)u{t,x,q,s) 
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fc + 7 V ' kJ Vq{t) 

Similarly, the maximizer for the ask part (for q ^ —Q) is: 



k \Vq-i{t)J 7 



and 



sup A" ((5") [u{t,x + s + 6°-,q - l,s) - u{t,x,q,s)] 

5» 



7 



fc + 7 



Aexp(-fc^''*)M(t, X, q, s) 



, 1 + 7) ' a;, g, s) 

fc + 7 V fc/ ^^(i) 

Henee, putting the terms altogether we get for \q\ < Q: 
dtu{t, X, q, s) + ^a'^d'^^u{t, x, q, s) 
+ sup A''(5'') \u{t, X — s + 6^,q + 1, s) — u{t,x, q, s) 

5'' 

+ sup A"(5") [u{t, X + s + 6'',q- l,s) - u(t, x, q, s) 
50. 



lVq{t) 7^0-^ 2 

— -u H — q u 



k Vq{t) 



-fA 
fc + 7 



7 u 



fc Vq{t) 



Vq{t) 



kja"^ 



= 



For q = —Q we have: 



dtu{t, X, q, s) + -a'^d1^u{t, x, q, s) 
+ supX^{6^) [u{t, X — s + 5^,q + 1, s) — u{t,x, q, s) 



-u + 



q u ■ 



7N ^ Vq+i{t) 



7 u 



kVq{t) 2 fc + 7 V fc/ Vq{t) 

v.it) - ^q\it) +a(i + lY^'^^^ Vq+,{t) 



kVq{t) 

Similarly, for g = Q we have: 



= 



dtu{t, x, q, s) + -(7^dl^u{t, x, q, s) 
+ sup X"'{5"') [u{t, X- s + S"-,q+l,s)- u{t, x, q, s)] 

Sa 



Dealing with ttie Inventory Risic 



27 



7 u 



k Vq{t) 



-k^f + ^''''-kT^[^+k) -^"^ 

Vq{t) - ^-^q\,{t)+A(l + ly^'^"^ Vq.,{t) 







Now, noticing that the terminal condition for Vq is consistent with the ter- 
minal condition for u, we get that u verifies (HJB) and this proves Proposition 
1. 

The positivity of the functions (^^g)|g|<Q was essential in the definition of 
u. Hence we need to prove that the solution to the above linear system of ordi- 
nary differential equations, namely v{t) = exp{—M{T — t)) x (1, . . . , 1)' (where 
M is given in Proposition 2), defines a family {Vq)\q\<cQ of positive functions. 

In fact, we are going to prove that: 

Vte [0,T].Vge {-Q,...,0}, > e-("«'-")(^-*) 

If this was not true then there would exist e > such that: 



min e"""^^ (vg{t) - e-(««'-'')('^-*)) + e{T - t) < 



^ -2v(T-t) 
te[0,T],\q\<Q^ 



But this minimum is achieved at some point {t*,q*) with t* <T and hence: 
d 



> e 
t=t* 



This gives: 

2^e-2''(^-**) (vq.{t*) - e-(«Q'-'')(^-**)) 

Hence: 

Now, if \q*\ < Q, this gives: 

aq*\.{t*) - v{vq^+i{t*) - 2vq.{t*) + V-i(i*)) 

-(77 + aQ2)g-(aQ=-r,)(T-r) > ^g2„(T-r) 

Thus: 

aq*' (^Vq, (f ) - e-(«0'-'')(^-**)) - v{Vq'+l{n - 2Vq. {t*) + Vq^-l{t*)) 
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All the terms on the left hand side are nonpositive by definition of {t*,q*) 
and this gives a contradiction. 

If g* Q, we have: 
Thus: 

-V{VQ-Iin ~ VQit*)) + + {vQin - e-(-Q^-"HT-*-)) > ,e'viT-n 

All the terms on the left hand side are nonpositive by definition of (t* ,q*) — 
{t*,Q) and this gives a contradiction. 

Similarly, if q* = — Q, we have: 

(aQ2 + r^)v^Q{t*) - v{v^Q+,{t*) - vgit*)) 
-(77 + ag2)e-("'3'-'')(^-**) > ee2''(^-**) 

All the terms on the left hand side are nonpositive by definition of {t* ,q*) — 
{t* , — Q) and this gives a contradiction. 

As a consequence, Vq(t) > e~("'3'^'')(^~*) > and this completes the proof 
of Proposition 2. 

Combining the above results, we see that u, as defined in Theorem 1, is a 
solution of (HJB). Then, we are going to use a verification argument to prove 
that u is the value function of the optimal control problem under consideration 
and prove subsequently that the optimal controls are as given in Theorem 1. 

Let us consider processes (z/'') and {v"-) G A. Let t G [0,T) and let us 
consider the following processes for r e [t,T]: 

dSy = adWr, S't' = s 

dXl'^^^" = (5, + v'^)dN^ ~ {Sr ~ yl)dNl, Xl'""'" = X 

dq','"'-' = dN^ - dN^, q't"^" = q 

where the point process has intensity (A^)r with — Ae^^'^^lq^_<^Q and 
where the point process N'' has intensity (A")i- with A° = Ae~^'^°lg, ^-nF^. 



These intensities are bounded since and u'^ are bounded from below. 
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Now, since u is smooth, let us write Ito's formula for u, between t and tn 
where tn = T A infjr > t, \Sr - s| > n or \N^ - N^\ > n or \N^ - N^\ > n} 
(n e N): 

u{tn, Xlr, qt'L', Sl'J = u{t, X, q, s) 

+ {u{t, X'^r - Sr + vl q^j" + 1, SY) - u{t, XYj\ qY-\ S'/)) X^r 
+ /*" (^x(t, xY-"" + SY + v'^, qY- - 1, SY) - u{t, X^^\ q^-. S^)) dM^ 



+ J'" {u{t, XY'" - SY + ^1 q^-" + 1, SY) - u{t, XYj\ q^- , S^)) dM\ 

where and M°- are the compensated processes associated respectively to 
and N°- for the intensity processes (A^)^ and (A"),-. 

Now, because each Vq is continuous and positive on a compact set, it has a 
positive lower bound and ?;g^(r)~^ is bounded along the trajectory, indepen- 
dently of the trajectory. Also, because and v"" are bounded from below, and 
because of the definition of tn, all the terms in the above stochastic integrals 
are bounded and, local martingales being in fact martingales, we have: 

E [u(t„, Xl'^:^ , ql^:^, Sl'^)] = u{t, X, q, s) 
+E (^dMr,XY''',qYj-,Stn + y ^^^(r, dr 

+ {u{t, X'^Y + SY + - 1, SY) - u(t, XYj\ qY-\ S^)) X^dr 

+ {u{t, X^J"" - SY + q^-" + 1, SY) - u{t, XYj\ qY-\ S^)) X\dT 



It 

Using the fact that u solves (HJB), we then have that 



E [u{tn, Xlr, iC-, S\:)\ < u{t, X, q, s) 

with equality when the controls are taken equal the maximizers of the hamil- 
tonians (these controls being in A because v is bounded and has a positive 
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lower bounded). 

Now, if we prove that 

lim E [u{U,, Xl^l\ 5*;^)] = E [«(T, X'^^-\ 4^)] 



we will have that for all controls in A: 



E [— cxp (— 7(X^^''' + qj'''''^ Sj.^))\ = E \u{T, X^'^'^, g^'"'', 5'^^)] < u{t, x, q, s) 
with equality for — d''*{t, qt-) and = S"-*{t, qt-)- Hence: 

sup E [— exp {—'y{X^^''' + g^^''^<S'^*'))] = u{t, x, q, s) 

= E [- exp (-7(X^"'^* + 4«'^*5^^)y 
and this will give the result. 

It remains to prove that 

lim E[u(^„,x*f_:^g*;«l^S*;^)] = E [u(T,x^-'^4«■^4^)] 

First, wc have, almost surely, that u{tn, X^'^z" ,qt'^- , S^'^) tends towards 
u(T, X!^^'^ , qj'-'^i S!j^''^). Then, in order to prove that the sequence is uniformly 
integrable we will bound it in L^. However, because of the uniform lower boimd 
on V already used early, it is sufficient to bound exp(— 7(X|'^'' + ql'''!^ Sl"'')) 

But, 



" " " Jt Jt Jt 



> - IKWooN^ - \\ut\\ooN^ + a J^'^ ql'^'^dWr 



Hence 



E [exp(-27(X|f_!- + qt'l" Sl':))] 



< E 



exp (27||z/» llooA/'^) exp {2j\\u^^N^) exp (-27<7 j[ " q^'^'^dW, 



< E [exp (67||i'^||ooiV^)] ' E [exp (67||i/^||ooiV^)] ' 



xE 



exp ( —6ja I q. 
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Now, since the intensity of each point process is bounded, the point pro- 
cesses have a Laplace transform and the first two terms of the product are 
finite (and independent of n). Concerning the third term, because is 
bounded by Q, we know (for instance applying Girsanov's theorem) that: 



E 



<E[exp(37V(t„-t)g2)]3 



<exp(7VQ2T) 



Hence, the sequence is bounded in L^, then uniformly integrable and we 
have: 

lim E[u(^„,x*f_:^g*;'^^5*:^)] =EKT,x^!■^9^!.'^5^^)] 

-E[u(r,X^"^",q^«'",5^^)] 

We have proved that u is the value function and that 6''* and 5"* are 
optimal controls. 

□ 

Proof of Theorem 2: 

Let us first consider the matrix M + 2ril. This matrix is a symmetric matrix 
and it is therefore diagonalizable. Its smallest eigenvalue A is characterized by: 



A= inf 



x'{M + 2r]I)x 



xm^Q+^\{o} x'x 
and the associated eigenvectors x ^ are characterized by: 

_ x'{M + 2r]I)x 
A — : 



It is straightforward to see that: 

Q Q-i 
x'{M + 2t]I)x = ^ aq^Xq^ + ^ X! ~ ^s)^ ^^^^ W-Q^ 

q=-Q q=-Q 

Hence, if x is an eigenvector of M + 2ril associated to A: 

^ ^ |a;|'(M + 2r;7)|x| 

A ; — r-j — ; 



1 



x'x 



Q 



q=-Q 



2|„ |2 



q=-Q 
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1 



a; ' a; 



Q Q-i 

^ aq^\xq\^ + 'q ^ {xq+i - Xqf + 'q\xQ\^ + 'q\x-Q\ 

q=-Q q=-Q 



X 



This proves that \x\ is also an eigenvector and that necessarily Xq+i and 
Xq are of the same sign {i.e. XqXq+i > 0). 

Now, let a; > be an eigenvector of M + 2r]I associated to A. 

If for some q with j^l < Q we have a;^ = then: 

= XXq = aq'^Xq - r]{Xq+l - 2Xq + Xg-l) = -T]{Xq+l + Xq-l) < 

Hence, because a; > 0, both Xq+i and Xq-i are equal to 0. By immediate 
induction a; = and this is a contradiction. 

Now, if a;Q = 0, then = Xxq = aQ^XQ-r]{-2xQ+XQ-i) = -r]XQ-i < 
and hence xq-i = 0. Then, by the preceding reasoning we obtain a contradic- 
tion. 

Similarly if x^q = 0, then = Xx^q = aQ'^x^q - ri{x^Q+i - 2x^q) = 
—r]X-Q+\ < and hence X-q+i = 0. Then, as above, we obtain a contradic- 
tion. 

This proves that any eigenvector a; > of M + 2r]I associated to A verifies 
in fact X > 0. 

Now, if the eigenvalue A was not simple, there would exist two eigenvectors 
X and y of M + 2r]I associated to A such that \x\'y = 0. Hence, y must have 
positive coordinates and negative coordinates and since yqVq+i > 0, we know 
that there must exist q such that yq = 0. However, this contradicts our pre- 
ceding point since \y\ > should also be an eigenvector of M + 2r]I associated 
to A and it cannot have therefore coordinates equal to 0. 

As a conclusion, the eigenspace of M -|- 2r]I associated to A is spanned by 

a vector /° > and we scaled its M^'^+^-norm to 1. 

Now, because M is a symmetric matrix, we can write v{0) = exp(— MT) x 
(l,...,l)'as: 



2Q 

=Y.exp{^X'T){g\ (1, . . . , l)')^;, Vg e {-Q, ...,Q} 

where A" < A^ < . . . < A^*? are the eigenvalues of M (in increasing order 
and repeated if necessary) and (5*), an associated orthonormal basis of eigen- 
vectors. Clearly, we can take = f^. Then, both and (/°, (1, . . . , 1)') are 
positive and hence different from zero. As a consequence: 



Dealing with ttie Inventory Risk 



33 



v,{0) ~T^+oo exp(-AOT)(/0, (1, . . . , 1)')/°, e {-Q, ...,Q} 

Then, using the expressions for the optimal quotes, we get: 

^ii„^s'.(0,„ = li„(i + 2) + ii.(|-) 

J,„^*-»(0,,) = il„(l + 2)+lM(|-) 

Turning to the eliaractcrization of /"^ stated in Theorem 2, we just need to 
write the Rayleigh ratio associated to the smallest eigenvalue of M + 277/: 

/°e argmin f'{M + 2r]I)f 

/eK2<3+Mi/|i2=i 



Equivalently: 

Q Q-i 

□ 



f G argmin ^ "^'/s' + ^ E (Z^+i - A)' + Vfg^ + Vf-Q^ 



Proof of Proposition 3: 

Let us first introduce H = {u e Lj^^{R)/x 1-^ xu{x) e L^(M) and u' € 



H equipped with the norm \\u\\h = y /r {oix'^u{xy + r]u'{x)'^) dx is an Hilbert 
space. 

Step 1: H C i^(M) with continuous injection. 
Let us consider u £ H and e > 0. 
We have: 

/ u{x)'^dx < -J / x'^u{x)'^dx < +00 

iM\[-e,e] e iK\[-e,e] 

Hence because u' G _L^(]R), we have u G H^(E. \ [— e,e]) with a constant Ce 
independent of u such that ||u||i/i(R\[_g^g]) < Ce||M||//. In particular u is con- 
tinuous on M*. 

Now, if e = 1, Va; G (0,1), = - j\'{t)dt and then |M(a;)| < 
\u{l)\ + VT^\\u'\\Lmo,i))- 
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Because the injection of _ff^((l,+oo)) in C([l,+cx))) is continuous, we 
know that there exists a constant C independent of u such that |w(l)| < 
C'||w||_f/i((i,+oo))- Hence, there exists a constant C such that |w(l)| < C'IIwHh 
and eventually a constant C" such that ||w||j;,oo((o,i)) < C"||u||ij. Similarly, we 
obtain ||w||loc((_i_o)) < C"\\u\\h- 

Combining the above inequalities we obtain a new constant K so that 

IKIk=(R) < K\\u\\h. 

□ 

A consequence of this first step is that H c H^{R) C C(M). 
Step 2: The injection H L^(K) is compact. 

Let us consider a sequence of functions in H with sup„ ||wn||// < +cxd. 

Because H C H^{R), Vm € N*, we can extract from (u„)n a sequence 
that converges in L'^{{—m,m)). Using then a diagonal extraction, there exists 
a subsequence of (■u„)„, still denoted (wn)n, and a function u G Lj^^iW) such 
that Un{x) — > u{x) for almost every x e M and u„ — )• u in the ^^^^(IR) sense. 

Now, by Fatou's lemma: 



/ a;^u(a;)^da; < liminf / x^Un{x)'^dx < 



SUPn \\Un\\H 

Hence, there exists a constant C such that Vm e N*: 



/ \u{x) — Un{x)\'^dx < [ \u{x) — Un{x)'\^dx-\--^ ( x'^\u{x)—Un{x)'\^ dx 

Jr J-m Js\[-m,m] 

/m Q 
\uix) - Un(x)fdx H 5- 



Hence limsup„^^ \u{x) - u„{x)\^dx < ^. 
Sending m to +oo we get: 

limsup / \u{x) — Unixy^ dx = Q 

n^oo Jr 

Hence converges towards u in the L^{R) sense. 

□ 

Now, we consider the equation —r]u"{x) + ax'^u{x) = f{x) for / e L^{M) 
and we define u = Lf the weak solution of this equation, i.e.: 

G H, (^ax^u{x)v{x) +r]u'{x)v'{x)) dx = / f{x)v{x)dx 
Jr Jr 
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Step 3: L : L^(R) L^(M) is a well defined linear operator, compact, 
positive and self-adjoint. 

For / G L^(R), V E H 1-^ J^f{x)v{x)dx is a continuous linear form on H 
because the injection H ^ L^(R) is continuous. Hence, by Lax-Milgram or 
Riesz's representation theorem, there exists a unique u G H weak solution of 
the above equation and L is a well defined linear operator. 

Now, = {f,Lf) < \\f\\mR)\\Lf\\L^R)- Hence, since the injection 

H ^ L^(R) is continuous, there exists a constant C such that < 
C||/l|L2(R)||L/||i/, which in turn gives < C||/||i2(R). Since the injec- 

tion H ^ L^(E) is compact, we obtain that L is a compact operator. 

L is a positive operator because {f,Lf) = > 0. 

Eventually, L is self-adjoint because V/, g e Z/^(R): 

(/, Lg) = [ {ax^Lf{x)Lg{x) + Tj{Lf )\x){Lg)'{x)) dx 
Jm. 

= [ {ax^Lg{x)Lf{x) + 7^{Lg)' {x){Lf)' {x)) dx = (<?, Lf) 
Jr 

a 

Now, using the spectral decomposition of L and classical results on Rayleigh 
ratios we know that the eigenfunctions / corresponding to the largest eigen- 
value A'^ of L satisfy: 

1 _ H/IU Mh 

A° ||/IU2(M) aeH\{0} MIl^k) 

Hence, our problem boils down to proving that the largest eigenvalue of L 
is simple and that g : x ^ exp is an eigenfunction corresponding 

to this eigenvalue (it is straightforward that g £ H). 

Step 4: Any positive eigenfunction corresponds to the largest eigenvalue of 

L. 

By definition of 11 • \\h, V/ € H, .'''{l"^ = JP^" . Hence, if / is an 

eigenfunction of L corresponding to the eigenvalue A*^, then |/| is also an 
eigenfunction of L corresponding to the eigenvalue A" . Now, if / is an eigen- 
function of L corresponding to an eigenvalue A ^ A°, (|/|, /) = 0. Therefore / 
cannot be positive. 

□ 

Step 5: g spans the eigenspace corresponding to the largest eigenvalue of L. 
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Differentiating g twice, we get g"{x) = —.^J^g{x) + ^x'^g{x). 

Hence —rjg"{x) + ax'^g(x) = y^arjg(x) and g is a positive cigcnfunction, nec- 
essarily associated to the eigenvalue A° that is therefore equal to ;^/=- 

Now, if wc look for an eigenfunction / G C°°(IR) OH- because any cigcn- 
function of L is in C°^(M) - we can look for / of the form / = gh. This 
gives: 

= -vf'ix) + ax''f{x) - V^fix) 
= -ri {g"{x)h{x) + 2g'{x)h'{x) + g{x)h"{x)) + ax^g{x)h{x) - ^Jotqg{x)h{x) 



Hence: 



= 2g'{x)h'{x) + g{x)h"{x) = -2x^g{x)h'{x) + g{x)h"{x) 



h"{x) = 2x^-h'{x) 



3K^, h' {x) = Ki exp ^x^ 
3Ki,K2, h{x) = Ktj^ exp (^^t^^ dt + K2 
3Ki,K2, fix)^Kig{x)J^ exp (^^t"^^ dt + K^gix) 



Now, 



9{x) r 

Jo 



exp ( .l—f] dt > exp ( -^.f^x'^ ] I exp [ .1 —t'' ] dt 



V J \ 2y T] J J. \\ T] 

Hence, for / to be in H, we must have Ki = 0. Thus, g spans the eigenspace 
corresponding to the largest eigenvalue of L and Proposition 3 is proved. 

□ 



